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Tunnel ionization belongs to the fundamental processes of atomic physics. The so-called two-step model, which 
describes the ionization as instantaneous tunneling at the electric field maximum and classical motion afterwards 
with zero exit momentum, is commonly employed to describe tunnel ionization in adiabatic regimes. In this 
contribution, we show by solving numerically the time-dependent Schrodinger equation in one dimension and 
employing a virtual detector at the tunnel exit that there is a nonvanishing positive time delay between the electric 
field maximum and the instant of ionization. Moreover, we find a nonzero exit momentum in the direction of the 
electric field. To extract proper tunneling times from asymptotic momentum distributions of ionized electrons, it 
is essential to incorporate the electron’s initial momentum in the direction of the external elechic field. 

PACS numbers: 03.65.Xp, 32.80.Fb 


In 1932, MacColl fT] studied the time which may he associ¬ 
ated with the process of a particle approaching from far away 
a potential barrier of a height larger than the particle’s energy 
and eventually tunneling through the barrier. Many efforts have 
been directed toward defining ||2]-0 and measuring ||3}0 tun¬ 
neling times. A related open problem is the question of how 
long it takes to ionize by tunneling through a binding poten¬ 
tial. Keller et al. have conducted experiments using the angular 
streaking technique aiming to measure tunneling times for ion¬ 
ization from a bound state, the so-called attoclock experiments 
CTOlfTTI . In the tunnel ionization case, a Coulomb-bound elec¬ 
tron is ionized by a strong electromagnetic field, and a potential 
barrier can be defined via the electron’s binding energy and 
the Coulomb potential bent by the electric field’s potential; see 
Fig-B Since the attoclock experiments have been performed, 
many renewed efforts have been directed toward defining a 
tunnel ionization time, because a consensus on a suitable the¬ 
oretical definition of tunneling time and the interpretation of 
experimental results is still lacking lHOHIIl- 

In the following, we study the time delay between the 
instant of ionization, i. e., when the electron exits the barrier, 
and the instant of electric field maximum fp. One may recon¬ 
struct the moment of ionization from the electron’s asymptotic 
momenrnm p{oo), which reads with the electron’s charge q 
and the time-dependent electric field E(t) (applying the dipole 
approximation and neglecting Coulomb corrections) 

p(oo) ^ p(to + ta) + q E(t')dt'. (1) 

*Jro+TA 

In attoclock experiments, electrons are ionized by elliptically 
polarized light, which makes the direction of the asymptotic 
momentum very sensitive to the ionization time (q + ta- In the 
extraction of ta from attoclock experiments, it is state of the 
art to treat the ionized electron classically, to take into account 
Coulomb corrections, and to assume that the electron’s initial 
momenrnm follows from some semiclassical theory naiisi. 
For a reliable reconstruction of the attoclock time ta, however, 
suitable initial conditions have to be identified as pointed out 
in Refs. ll20ll2Tll . In particular, assumptions about the initial 
momentum bias the reconstructed value for ta- 


The popular two-step model of tunnel ionization assumes a 
maximal ionization rate at the instant of maximal electric field 
strength, i. e., ta = 0, and that free electrons have zero initial 
momentum, i. e., p{to + ta) = 0. Within the two-step model, 
the electron’s asymptotic momentum follows by solving the 
classical equations of motion for the electron’s motion in the 
combined electromagnetic field of the binding potential and 
the ionizing external light pulse. This model, however, cannot 
be used as a benchmark for experiments, because a possible 
match or mismatch of experimental data and the theoretical 
prediction by the two-step model may be explained by various 
pairs of nonzero delay ta ll22ll and nonzero initial momentum 
p{to + Ta), which may a have nonzero component parallel to 
the electric field direction ll^l24ll . 

Neither the delay ta nor the initial momentum p(to + ta) are 
directly accessible by experiments, and it is also challenging 
to calculate them analytically. Therefore, we employ ab initio 
quantum calculations and a virtual detector Il25l l26ll at the tun¬ 
nel exit. The virtual detector technique allows us to determine 
directly the electron’s time of arrival at the tunnel exit as well 
as its exit momentum. 

We analyze theoretically an initially bound electron ionized 
by an electric field pulse. To determine the time delay ta and 
the exit momentum, we solve the time-dependent Schrodinger 
equation and place a virtual detector at the tunneling exit. The 



FIG. 1 : The one-dimensional Coulomb potential, which is bent by an 
additional electric field, defines the effective potential V(x) and the 
tunneling barrier between the points and x^xit for an electron, which 
is given by the wave function ij/oix), with the ground state energy -Ip. 
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virtual detector is realized by calculating the probability current 
at the exit. The most probable time delay ta is determined 
by comparing the instant of maximum probability current at 
the potential barrier exit and the instant of maximum electric 
field strength. The exit momentum is determined by the space- 
resolved momentum distribution at the tunnel exit at the instant 
of ionization. By separating the wave function into a tunneled 
part and a bound part after the interaction with the laser pulse, 
we can calculate the momentum distribution of the tunneled 
electron, from which one can determine the most probable 
asymptotic momentum. 

Considered system In experiments, a Coulomb-bound 
electron is usually excited by a laser pulse with a wavelength 
much bigger than the atomic dimensions such that the laser 
pulse is nearly homogeneous over the size of the atom. Further¬ 
more, relativistic effects and effects due to the magnetic field 
component set in only for tunneling from highly charged ions 
EU- Thus, we will apply the electric dipole approximation. 
The laser pulse is modeled by a time-varying homogeneous 
electric field E{t) - Eq exp(-m^(f - to)^l2), where Iq denotes 
the instant of the maximum field strength Eq and te = V2/m 
is the time scale of the rise and decay of the electric field. 
The linear polarization of the electric field and neglecting the 
magnetic field component render the motion of the electron 
quasi-one-dimensional allowing us to investigate general fea¬ 
tures of tunneling times in a one-dimensional scenario. Fur¬ 
thermore, tunneling in the fully three-dimensional Coulomb 
problem can be described by an effective one-dimensional tun¬ 
neling barrier via introducing parabolic coordinates ll28l . Thus, 
we restrict ourselves to one-dimensional systems and consider 
an electron bound to the soft-core potential -Z/ 

12914321 to model the essential features of an electron in a three- 
dimensional Coulomb potential. Here, Z is the atomic number, 
and the softening parameter a(Z) = 2/Z^ is chosen such that the 
ground state energy of the soft-core potential is -/p = -Z^/2, 
which equals the ground state energy of the Coulomb potential 


X = Xin - = Xin + 2Ax - JC = JCjn + 4Ax 

x = JCin + AjC - - = Xin + 3Ax - ^ = ^exit 



(t - to)Z^ (a.u.) 


FIG. 2: The probability current j{x, t) as a function of time at different 
positions x between the barrier entry x = Xjn and the barrier exit Xe^it 
separated by Ax = (Xe^it -Xi„)/5 for the parameters Eq/Z^ = 0.048 and 
7 = 0.25. The shadowed area represents the time between the instant 
of maximum electric field strength tp and the instant of maximum 
probability at Xexit- 


Ea. Thus, the Schrodinger equation (in atomic units) 

Z 


dtp 

= H{t)tp = 
at 


1 

2 dx'^ 




■a{Z) 


- E(t)x 


iP ( 2 ) 


with the Hamiltonian H(t) and the effective potential V(x, t) — 
-Zj xjx^ + a{Z) - E(t)x will be solved numerically ll34l . The 
so-called Keldysh parameter y - at -sJlLplEo fia characterizes 
the ionization process as dominated by mnneling for y « 1 
and by multiphoton ionization for y » 1. Thus, simulation 
parameters will be set such that y < 1 in the following. 

Time delay ta The time delay ta is based on the 
time-dependent ionization rate. In our one-dimensional 
model the probability current j{x,t) - {tp{x,tydxip{x,t) - 
ip{x, t)dxip{x, f)*)/(2i) represents the average net number of 
electrons passing a given point at a specific time. Thus, we 
determine the ionization rate via the probability current at the 
exit Xexit as a function of time, where Xexit is defined by the 
maximum electric field strength Eq via y(xexit, to) - -Ip. Mon¬ 
itoring the probability current at a fixed position is justified, 
because the tunnel probability is maximal for E{t) - Eq and it 
is exponentially suppressed for lower electric fields. Further¬ 
more, E(t) - Eo{l - + (9(Af"^/Tg)) for At - t - to with 

|Af| < Te = 2 VTp/Cyfio)- Thus the tunneling barrier does not 
change substantially if times close to to are considered. To 
further ensure that the calculated exit Xexit is where the electron 
exits the barrier, we place a virtual detector at different points 
between Xin and Xexit and calculate j{x,t) at each point as a 
function of time. As shown in Fig. the probability current 
has a positive peak as well as a negative one for x < Xexit, in¬ 
dicating the tunneling and reflection dynamics, respectively. 
This tunneling and reflection dynamics corresponds to under- 
the-barrier dynamics. As reflection is absent for x > Xexu, the 
particle leaves the barrier at Xexit- 

The most probable time delay ta is calculated by subtracting 
the instant of maximum field strength from the instant of the 
maximum current, i. e., ta = arg max /(Xexit, t) - arg max £(?), 
which yields the positive time delay shown in Fig.[^ Note that, 
for the parameters used in Fig.[^ ta < te, and thus the electric 
field remains almost constant for times |Af| ^ ta, justifying our 
choice for Xexit- 

The origin of the time delay ta can be understood by consid¬ 
ering the time-energy uncertainty principle and following its in¬ 
terpretation given by Mandelstam and Tamm lfT4l 13514371 . As a 
consequence of the time-energy uncertainty principle, the time, 
which a wave function of a system with a time-independent 
Hamiltonian H needs to change significantly, is bounded from 

below by l/(2cr^), where cr^ = x^{ip\H'^\ip) - {tp\H\>p)'^ . As 
outlined in the previous section, the Hamiltonian in Q can 
be considered as time-independent for times |Af| < te. This 
allows us to define the Mandelstam-Tamm time 


1 / - - 
tmt = 2 (<*A(fo)|H(fo)^|tA(fo)> - {tp{to)\H{to)\p{to)) j 

(3) 

which is indeed a lower bound to the time delay ta as indicated 
in Fig. The time delay ta is close to its lower bound tmt> 
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• Tk,y = 0.25 ♦ tmt, r = 0.25 a ri,y = 0.25 — Ts„b 

o Tk,y = 0.35 o Tmt. y = 0.35 A Ti,y = 0.35 



FIG. 3: The time delay Ta, the Mandelstam-Tamm time Tmt, and the 
times T 2 and Tsub plotted for different electric field strengths £0 and 
for different Keldysh parameters y. Definitions of the various times 
are given in the main text. 

which indicates that the delay ta is a consequence of the wave 
function’s inertia, i. e., its inability to adopt instantaneously to 
the field. The time delay ta decreases as Eo increases at fixed 
Keldysh parameter j and matches approximately tmt when the 
regime of over-the-barrier ionization is approached, which is 
for EqIZ^ ss 0.06 a.u. 

The observed decrease of the delay ta with growing electric 
field strength (but constant y) is consistent with calculations 
for the case of a sudden turn-on of the electric field, which 
show that the time for the wave function to adopt to the elec¬ 
tric field is proportional to the Keldysh time tk = sflLplEo 
01 EH, although the functional dependence on Eq is different 
for continuously changing electric fields. For a fixed maximum 
electric field strength Eq but increasing y, the time delay ta in¬ 
creases. As y increases, the pulse duration decreases, granting 
the wave function less time ~ te to evolve and to adopt to the 
time-varying Hamiltonian during the rise of the electric field. 
Thus, the wave function has less time to develop the necessary 
components for tunneling and thus needs more time to reach 
the maximum ionization rate. 

The times ta and tmt are intractable by analytical methods. 
As we will demonstrate in the following, however, one may 
employ the Wigner time as an approximation for ta. Compar¬ 
ing the quantum mechanical Wigner trajectory fwC^c) caEii 
to the trajectory t^x) of a classical particle, which travels in¬ 
stantaneously from Xjn to Xexit and then moves according to the 
classical equations of motion, we can calculate the time delay 

^sub “ ^w(.^exlt) C(.^exlt) . (4) 

where the Wigner and the classical trajectories fw(.*) and fc(x) 
are determined such that both coincide at the entry point Xin. 

The delay ta describes the time interval that the electron 
spends under the barrier after the field maximum. Interpreting 
Ta as a tunneling time IIT^ is justified only if the classical for¬ 
bidden region is entered at the field maximum. The Wigner 
delay 0. however, represents the total time spent under the 



/ (lOi'^W/cm^) 

FIG. 4: Wigner time Ts^b, which is an estimate for r^, for the three- 
dimensional Coulomb potential including Stark shift and multielectron 
effects as a function of the electromagnetic field’s peak intensity I 
for ionization of helium atoms by linearly polarized laser fields. The 
applied wavelength is T = 800 nm and the Keldysh parameter varies 
from 0.6 to 0.2. 

barrier, since the relation fc(.*in) = fc(jCexit) holds in the classical 
model. Also, the time delay ta is often understood as the time 
under the barrier. Our numerical results show that Tsub is indeed 
close to the time delay ta, as shown in Fig. where Tsub is cal¬ 
culated for the one-dimensional Coulomb potential. Although 
Tsub and Ta agree well, it is open if ta should be interpreted 
as time spent under the barrier, because the interpretation as 
under-the-barrier time is based on the assumption that tunnel¬ 
ing starts at the instant of the maximal electric field, but the 
wave function could penetrate the tunneling barrier earlier. 

For the one-dimensional Coulomb potential and at field in¬ 
tensities approaching the over-the-barrier threshold, the Wigner 
time is given by Tjub ~ 14.29 yjl - \6 EqIZ^IZ^, and by Tsub ~ 
9.0 yjl - 9.5£’o/Z^/Z^ for the three-dimensional case. Further¬ 
more, the Wigner formalism can be extended lIMl to include 
the Stark shift and also multielectron effects. The resulting 
Wigner time Tjub for helium is shown exemplarily in Fig. as 
a function of the laser’s peak intensity. 

Exit momentum In the following, we determine the quan¬ 
tum mechanical exit momentum pQ in the direction of the 
electric field at the tunnel exit Xexit at the instant of ionization 
f = fo -H Ta by two different methods. In the first method, we cal¬ 
culate a space-resolved momentum distribution around x = Xexu 
by multiplying the wave function fi{x, to + ta) by a Gaussian 
window function with mean Xexit and width 6x - (Xexit - .*in)/20 
and calculating its Fourier transform i^exit(fo + ta)- The most 
probable exit momentum po can be inferred by the momentum 
where It^exitC^o + ta)P is maximal. In the second method, the 
most probable initial momentum po is determined from the 
probability current y(Xexit, fo + ta) at the exit at the instant of 
the maximum probability current. Following Refs. ESIESI, 
the local velocity of the wave function’s probability flow at 
Xexit equals v = j'ixexn, to + TA)/\fiXexiu to + ta)P- Both meth¬ 
ods yield very similar results for the moment po as shown in 
Fig. 0 The initial momentum po is almost independent of 
the parameter y and depends only weakly on the electric field 
strength Eq. 

As Po does not depend on the parameters of the external 
electric field it must result from the initial quantum state, i. e., 
the ground state of the binding potential. In fact, the ground 
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FIG. 5: The momentum po at the time to + Ta at the tunnel exit for 
different elech'ic field strengths Eq and different Keldysh parameters 
7 as determined by two different methods. Method 1 is based on the 
space-resolved momentum distribution, while method 2 utilizes the 
velocity of the probability flow; see the main text for details on both 
methods. 


state of the employed soft-core potential has in momentum 
space a width of about 0.38 x Z, which is of the same order 
as po; see Fig.|^ The fact that po scales with the width of the 
ground state’s momentum distribution and not with its mean 
(which is zero) may be interpreted as momentum components, 
which propagate into the ionization direction, being ionized 
preferably. 

Implications The above considerations are relevant for 
every high-precision laser-induced tunneling experiment. Here 
we discuss as an example attoclock measurements due to their 
currently high attention. Attoclock experiments aim to deter¬ 
mine the time delay t\ between the instant of the electric field 
maximum and the instant of tunneling, which is not directly 
accessible experimentally. Instead, the asymptotic momentum 
of the tunneled electron is measured. As it depends on the exit 
momentum and on the moment at which the electron starts to 
propagate freely in the field, one can infer t\ from the elec¬ 
tron’s asymptotic momentum provided that the exit momentum 
is known. The delay ta is commonly reconstructed by assum¬ 
ing zero initial momentum in the electric field direction. Our 
numerical simulations and the virtual detector approach, how¬ 
ever, indicated a nonzero initial momenrnm in the direction 
of the electric field. How does the zero-initial-momentum as¬ 
sumption affect the reconstruction of ta from the asymptotic 
momentum? 

To answer this, we determine the final momentum of the 
tunneled electron by propagating the wave function till 
some final time tf such that t{ - to » l/« and separate the 
tunneled part ij/beeitf) from the bound part of the wave function 
by projecting out all bound eigenstates of H in Eq. (|^ for 
E{t) - 0 from From the resulting probability density 

in momentum space representation lAfree(h) the most probable 
momentum pf q can be inferred by the position where |(^free(tf)l^ 
is maximal. Using the Newton equations of motion for an 
electron in the effective potential V(x, f), we can calculate at 
which instant of time the electron must exit the barrier at ;iCexit 
with zero initial momentum in the electric field direction such 
that its asymptotic momentum equals pf_q. The result is the 


time delay T 2 , which is also shown in Fig.j^and that does not 
coincide with the delay ta- The delay T 2 is close to zero or 
even negative depending on the electric field strength as found 
in Refs. Il20l - l22l 12411271 . Consequently, the delay ta and the 
instant of tunneling cannot be determined on the basis of the 
standard assumption of zero initial momentum in the electric 
field direction. The exit momentum has to be included. 

Similarly to the delay T 2 , also the time delays determined 
from measured asymptotic momenta are close to zero (of the 
order of experimental uncertainties) IfTOllTTl . As our numer¬ 
ical calculations indicate that the reconstruction of the delay 
Ta from asymptotic momenta is sensitive to the electron’s exit 
momentum, we argue that a possible nonzero exit momentum 
into the electric field direction has to be included for a reli¬ 
able reconstruction of ta in attoclock experiments, not only a 
nonzero exit momentum in the direction perpendicular to the 
electric field as recently proposed EOl [39l . 

Conclusions We reexamined tunneling times in strong 
field ionization by an ab initio solution of the time-dependent 
Schrodinger equation. Our calculations show that there is a de¬ 
lay Ta between the maximum of the ionization rate and the maxi¬ 
mum of the electric field strength and a nonzero exit momentum 
po in the electric field direction. The delay ta can be explained 
as the response time needed by the wave function to react to the 
change of the driving electric field. The initial momentum can 
be estimated from the width of the ground state’s momentum 
distribution, which is Z for the three-dimensional Coulomb po¬ 
tential. The time may be estimated by the Wigner time Tjub- 
Note that our results differ from the vanishing tunneling delay as 
obtained by the recently introduced analytical R-matrix method 
ED, which is based on the assumption that this quasiclassical 
model provides a good description of the quantum-mechanical 
tunneling dynamics also in the vicinity of the tunneling barrier. 

The authors thank Karen Z. Hatsgortsyan, Michael Klaiber, 
Oleg Skoromnik, and Anton Wollert for valuable discussions. 
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